Ripples in a string coupled to Glauber spins 
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Each oscillator in a linear chain (a string) interacts with a local Ising spin in contact with a 
thermal bath. These spins evolve according to Glauber dynamics. Below a critical temperature, a 
rippled state in the string is accompanied by a nonzero spin polarization. The system is shown to 
form ripples in the string which, for slow spin relaxation, vibrates rapidly about quasi-stationary 
states described as snapshots of a coarse-grained stroboscopic map. For moderate observation times, 
ripples are observed irrespective of the final thermodynamically stable state (rippled or not). 
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Mechanical systems coupled to spins are used to de- 
scribe structural phase transitions. Examples include the 
collective Jahn- Teller effect structural phase transi- 
tions with a scalar order parameter exhibiting a central 
peak in the dynamic response function Q and critical- 
ity in martensites and externally driven models jl]. In 
many of these models, the mechanical system provides a 
long range interaction bctvifeen the spins that produces a 
phase transition in which the spin polarization ceases to 
be zero below a critical temperature. While most studies 
consider the effective spin system obtained after eliminat- 
ing the mechanical degrees of freedom, it is interesting to 
focus instead on the effect of the phase transition on the 
mechanical system. In this paper, we consider mechani- 
cal systems coupled to Ising spins that undergo Glauber 
dynamics [J] in contact with a thermal bath (a single har- 
monic oscillator connected to Ising spins in the simplest 
case [^). There is a phase transition at a critical temper- 
ature below which the spin polarization is nonzero and 
ripples appear in the mechanical system. These thermo- 
dynamically stable ripples are inhomogeneous stationary 
states of the mechanical system, which are quite sim- 
ple below the critical temperature. On the other hand, 
there are long-lived dynamical ripples with a wide va- 
riety of shapes at any temperature provided the period 
of mechanical vibrations is short compared to the spin 
relaxation time. In this limit, the spins are frozen dur- 
ing long time intervals between spin flips and they fix a 
quasi-stationary state about which the mechanical sys- 
tem oscillates. Observations of the system may consist 
of time averages over intervals sufRciently long to include 
many oscillation periods but short compared to the in- 
tervals between spin flips. Then these observations will 
sample a coarse-grained stroboscopic map consisting of 
successive quasi-stationary states that show ripples. Af- 
ter a much longer time during which sufficiently many 
spin flips have occurred and due to the dissipation intro- 
duced by the Glauber spin dynamics, the ripples eventu- 
ally evolve to the simple version obtained from the equi- 



librium thermodynamics of the spin-mechanical system. 

These considerations may apply to the evolution of 
ripples in suspended graphene sheets. Ripples are on- 
dulations of the sheet with characteristic amplitudes and 
wave lengths that, according to experiments, do not have 
a preferred direction Time resolved ripples and de- 
fects in graphene sheets can be observed using aberration 
corrected TEMs that collect data every other second, a 
time much longer than microscopic times such as the one 
it takes a sound wave to cross one lattice constant R. As 
a direct generalization of theories of defect motion in pla- 
nar graphene Q, atom motion in a suspended graphene 
sheet may be described by the von Karman equations 
discretized on a hexagonal lattice Coupling the ver- 
tical motion of graphene atoms with an Ising spin located 
at the same lattice point may account for a spontaneous 
trend of the sheet to bend upwards or downwards. Spin 
dynamics adds dissipation to the von Karman equations 
and thus the spin relaxation time should be much longer 
than microscopic mechanical times. Experimental obser- 
vations are taken over long time intervals and therefore 
should correspond to different takes of a coarse-grained 
stroboscopic map similar to that described in this paper. 

Model and rippling phase transition. Our mechanical 
system is a chain of oscillators with next-neighbor inter- 
action which becomes a string in the continuum limit: 
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uj^^i. The jth oscillator is coupled 



linearly to an Ising spin aj — ±1. The spins are in 
contact with a thermal bath at temperature T and flip 
stochastically following Glauber's dynamics [1] at tem- 
perature T. Unlike the case of a regular mass-springs 
chain, here each triplet is biased against being straight, 
with Uj+i — Uj = Uj — Uj-i. The applied force (whose 
sign flips at random) makes the "preferred" state at any 
instant a wedge shape. This translates the loose idea 
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that the three carbon links each atom shares to build the 
graphene sheet do not want to be in a plane because of the 
fourth "free" link (which may push the atoms up or down 
the horizontal planar configuration). The free chemical 
bonds of the carbon lattice in the graphene sheet may 
be assimilated to our spins, which interact with phonons 
modeled by the oscillators. Of course, in order to have a 
more realistic model of a graphene sheet, the structure of 
its 2D lattice should be taken into account. Nevertheless, 
we hope that this simple model will be able to capture 
the main physical mechanism involved in the rippling of 
graphene sheets. 

Thus at any time t, the system may experience a tran- 
sition from (u, p, a) to (u, p, Rja) at a rate given by Q 
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where R^a is the configuration obtained from a by flip- 
ping the j-th spin and is the Boltzmann constant. 
The parameter a gives the characteristic attempt rate 
for the transitions in the Ising system. Individual spins 
experience a mutual long-range interaction through their 
coupling to the string. This long-range interaction causes 
a phase transition of the spin system: the spins have 
non-zero polarization for T < Tc whose counterpart 
for the string is the formation of ripples. To see this, 
we find the following effective potential by integrat- 
ing e"^/*^'^^^-' over the spin configurations [5|: VeS = 

Ef=o{"^w'("j+i - %)V2 - fcsrincosh[/M,/(A;BT)]}. 
The extrema of this potential satisfy 
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which, linearized about the trivial solution Uj — (hori- 
zontal string) yields 
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for j — 1, . . . N with boundary conditions Uq — Un+i — 
0. At T = Tc, (HI should have a nonzero solution 
Uj = e'^''. Inserting this in (g]), we get p / {mto'^kBTc) = 
4sin^(A;/2). A linear combination of e'-''^ and e"*-''^ that 
satisfies ?7o = is Uj ~ sin(jfc). The other boundary con- 
dition Un+1 — yields fc„ = m:/{N -|- 1), n = 1, . . . ,iV. 
The largest possible critical temperature corresponds to 
n = 1 and therefore 
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The critical temperature Tc remains finite as — > 00 
if we define a finite such that uj — ujoKn. At T = 
Tc = /^/ {muiQkB), there is a phase transition so that the 
trivial solution of ([3]) is linearly stable for T > Tc and 
unstable otherwise. We introduce the nondiniensional 



variables u* 



fuj/{kBTc) = muj^uj/f, set 9 = T/Tc 



and omit the asterisks so as not to clutter the formulas. 
Below the critical temperature, stable non-uniform states 
corresponding to ripples in the string appear. The first 
such state is 
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as the reduced temperature 6 = T /Tc tends to 1 from 
below. This result follows from bifurcation theory for 
the macroscopic equations for the averages {uj) and {uj) 
p^ . At lower critical temperatures Qn = sin^ ki/ sit? kn, 
n > 1, other non- uniform states with n ~ I interior 
nodes bifurcate from the trivial solution. The non- 
uniform solutions have a nonzero spin polarization (aj) ~ 
tanh(/Mj/(fcBT„)) and therefore the critical temperature 
Tc is associated with cooperative Jahn- Teller phase tran- 
sitions, in which coupling to phonons (the string) breaks 
the symmetry of a doubly degenerated electronic state 
(the spins) Numerical simulations of the spin-string 
system confirm this. Below Tc, stable string configura- 
tions are stationary, nonuniform and exhibit ripples. To 
test the bifurcation theory, we have performed stochas- 
tic simulations at temperatures 6 > 1 and 9 = 0.9 for 
S — 0.1 and N — 10*. In Figure [1] we show how an 
initially flat string at rest evolves to a state close to 
The initial conditions are random spins, such that the 
average spin polarization has a sinusoidal shape, and a 
horizontal, zero-velocity string profile. A qualitatively 
analogous behavior is found for the majority of initial 
spin configurations, but for some of them the unstable 
flat string configuration is stabilized ^Td\ . This is a stabi- 
lization of the thermodynamically unstable state akin to 
the one previously found for a single oscillator coupled 
to Glauber spins Q. For lower temperatures a similar 
stationary state without internal nodes is stable whereas 
the stationary states with n — 1 internal nodes that bi- 
furcate from the flat string configuration at temperatures 
T„ are unstable. 

Slow spin relaxation. Additional insight can be ob- 
tained in the limit S = a/ojo <C 1, in which the spin 
flip rate is slow compared to the characteristic string fre- 
quency. In the following, we will use a dimensionless time 
t* = u!ot, and omit the asterisks as before. Then the Uj's 
obey the equation of motion 
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for j = 1, . . . , N with boundary conditions uq = UAf+i — 
{il — d?u/dt^). The spins aj are stochastic variables 
which flip at a rate W"j((t|u, p) = 5{1 — /3jO'j)/2, jij — 
tanh(uj/6'), 9 = T/Tc instead of Let us consider a 
trajectory of the system, for given initial states of the 
string and spins. Since the spin flip rate is very small, 
the spins are frozen at fixed values during time intervals 
that are long compared to the longest oscillation period 
of the string. During the time interval before the spin 
flip occurs, we may split the solution of ([7]) in a quasi- 
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FIG. 1: Stable ripple state without internal nodes for 6 — 
0.9, S = a/u!o ~ 0.1 and A*' = 10^. For each trajectory, 
the initial configuration consists of a fiat string at rest and 
random spins, such that the average spin polarization has a 
sinusoidal shape. Averages over 100 trajectories and spatial 
averages over 100 oscillators centered at a given one have been 
performed in order to ensure good averages. See also movie 
in the Supplementary Material. 



stationary and a time-dependent part according to [lo| : 
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The string profiles ([SJ represent vibrations of the 
string about the quasi-stationary configuration Q whose 
longest period is 27r (ili = 1 is the lowest frequency). 
Now let the first spin that fiips after t ^ he Cj^ , 
which changes sign at time ti. Immediately after ti, 
the right hand side (RHS) of ([7]) should be replaced by 
aj-26j,jaj^Q{t-ti), where e(a;) = 1 for x > 0, Q{x) = 
for a; < is the unit step function. The changes in Uj 
and Vj{t) due to the spin fiip at ti are 
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respectively, for t > ti. Successive spin fiips produce 
changes similar to (IHl) and in the quasi-stationary 
and time-dependent parts oiuj{t), respectively, at times 



i2, h,...wit\iti-ti_i =0((iV5)-i). Time averages over 
sufficiently long time intervals that are short compared 
to {N6)^^ eliminate Vj{t). Thus successive snapshots 
of averaged string profiles coincide with updated quasi- 
stationary profiles. The latter constitute a coarse- 
grained stroboscopic map showing how the ripples in the 
string evolve to their final stable configurations: the hor- 
izontal string for > 1 or a simple parabolic-like profile 
above or below the horizontal string for < 1. Figure [2] 
depicts snapshots of the coarse-grained stroboscopic map 
for an initially flat string at rest with (a) a spin configu- 
ration exhibiting seven domains and (b) completely ran- 
dom spins for a temperature 9 — 0.1, below the critical 
one. For moderate time intervals stable ripples are ob- 
served whereas the stationary configuration of the string 
without internal nodes is reached at extremely long time 
intervals. Rippling behavior is also observed for above 
critical temperatures, 9 > 1, but the string eventually 
approaches the flat configuration. In time-resolved ex- 
periments such as those with suspended graphene sheets 
[7|, data are taken in long time intervals (1 second), which 
we use as an estimate for the spin fiip attempt rate a. 
Typical times w^^^ are 1 ps, therefore 5 ~ 10~^^ and the 
number of spins per linear dimension ~ 10* for square 
1 micron samples. Thus, \/{N5) ~ 10* (much larger 
than the value considered in Fig. [2]) and ripple states 
corresponding to snapshots of the coarse-grained stro- 
boscopic map are observed. The "true" thermodynam- 
ically stable state would only be reached for extremely 
long time intervals, much longer than the total observa- 
tion time in an experiment. On a physical basis, one 
may expect that the characteristic time associated to the 
spin fiips in graphene be larger than the data-collecting 
time (1 s), so our estimate for the time between spin 
fiips is actually a lower bound to the actual value. The 
movie in the Supplementary Material illustrates how the 
string vibrates rapidly about the quasi-stationary con- 
figurations corresponding to successive snapshots of the 
coarse-grained stroboscopic map. Given the large sepa- 
ration between microscopic times, data collection times 
and duration of a given experiment, it is important to 
remark that ripples are observed for all current time- 
resolved experiments no matter what the temperature and 
the thermodynamically stable state are. Thus ripples are 
inherently dynamical and explanations based on thermo- 
dynamically stable states do not capture the essence of 
rippling. 

Continuum limit and fast spin relaxation. Further 
analysis of string ripples can be done in the continuum 
limit N —i' oo. From ([T]) and ([2]), we obtain equations 
for averages of Uj and aj. If we split the variables 
Uj — Uj -\- A.Uj, where Uj = {uj), set qj — (cj) and 
ignore the fluctuations Amj in the limit as — > oo, we 
get the following nondimensional macroscopic equations: 
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FIG. 2: Four snapshots of the coarse-grained stroboscopic 
map for 9 = 0.1, S = 10"" and N = 100 up to t = 10^ 
The initial configuration consists of a flat string at rest and 
the spins are (a) initiaUy distributed in seven spin domains of 
alternating sign or (b) randomly distributed. The times cor- 
responding to the spin flips which change the quasi-stationary 
string proflle are indicated. Ripples with several domains 
persist for long times whereas the stable configuration cor- 
responding to a string without internal nodes is reached in 
a much longer time larger than 10*. See also movie in the 
Supplementary Material illustrating how the string oscillates 
about quasi-stationary configurations given by the coarse- 
grained stroboscopic map. 



for j — 1,...,N. We now set Uj{t) = u{x,t) with 
X = j/K^ and take the continuum limit. Then (|16|) - 
([T7| become 
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to be solved with the boundary conditions u(0, t) = 
u{7T,t) = 0. In the limit (5 3> 1 (fast relaxation of the 
spins compared to the string time scale) , we can approx- 
imate the second equation in (|T8| by q « ta.nh{u/6) + 



[S9 cosh^{u/ 9)] ^du/dt and insert this in the first equa- 
tion. The result is 
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whose stationary solutions are ^ just below = 1. The 
small damping term in (|19p stabilizes the trivial solution 
above the critical temperature and the stationary ripple 
solutions below it. The opposite limit of (5 <C 1 has al- 
ready been studied using the coarse-grained stroboscopic 
map. 

Conclusions. We have shown that stable ripples ap- 
pear in a ID string when each oscillator is coupled to an 
Ising spin and the latter are in contact with a thermal 
bath at temperature T. Below a critical temperature, 
the thermodynamically stable string profile is not flat, 
but nonuniform without internal nodes. In spite of the 
simplicity of the thermodynamically stable state, more 
complex ripples appear when the spin flip rate is much 
smaller than the oscillator period. Although strictly 
speaking these ripples are evolving in time, they are very 
long-lived metastable states. The ripples are snapshots 
of a coarse-grained stroboscopic map depicting the aver- 
age of the rapid string motion over long time intervals. 
Whether the final thermodynamically stable is the fiat 
or bent string, ripples should be observed on reasonable 
time intervals at any temperature. 

The system considered here is far from being a realis- 
tic model of a graphene sheet. However, the free chem- 
ical bonds of the carbon lattice in the latter may be as- 
similated to our spins, which interact with the phonons 
modeled by the oscillators. Thus, 2D ripples analogous 
to the ones found here should appear. This opens the 
door to understanding the characteristic rippling shown 
by graphene sheets at any temperature as an inherently 
dynamical phenomenon, whose physical mechanism con- 
sists of the interaction between free bonds and phonons 
with widely separated time scales. 
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